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Abstract 

We present a re-interpretation of the recent ATLAS limits on supersymmetry in channels with jets (with 
and without b-tags) and missing energy, in the context of Ught third family squarks, while the first two 
squark families ai^e inaccessible at the 7 TeV run of the Large Hadron Collider (LHC). In contrast to inter- 
pretations in terms of the high-scale based constrained minimal supersymmetric standard model (CMSSM), 
we primarily use the low-scale parametrisation of the phenomenological MSSM (pMSSM), and translate 
the limits in terms of physical masses of the third family squarks. Side by side, we also investigate the lim- 
its in terms of high-scale scalar non-universality, both with and without low-mass sleptons. Our conclusion 



O 
O; 

\ is that the limits based on zero-lepton channels are not altered by the mass-scale of sleptons, and can be 



considered more or less model-independent. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) at CERN is now running and the two experiments ATLAS 
and CMS are hard at work in constraining theories beyond standard model (BSM). Supersymme- 
try (SUSY) yj, |2D has long been one of the most popular BSM models, due to its ability to solve 
the hierarchy problem and to provide a dark matter candidate in its R-parity conserving versions. 
The minimal supersymmetric standard model (MSSM) extends the nineteen standard model pa- 
rameters to over a hundred. However, the proposed mechanisms of SUSY breaking, together with 
the requirements of suppresing flavour-changing neutral currents (FCNC) and obtaining the cor- 
rect electroweak breaking scale often suggest a common origin of some of the parameters at high 
scale. This results in highly constrained models like the minimal supergravity model (mSUGRA) 
1 3, 4] or the very similar constrained MSSM (cMSSM) [5], gauge mediated SUSY |^ and others. 
Alternately, the new parameters can be constrained by demanding only that all observed experi- 
mental constraints are obeyed, thus arranging SUSY breaking mass parameters in such a way as 
to suppress flavour-changing neutral current (FCNC) processes. This leads to a phenomenological 
MSSM (pMSSM) |9] which has nineteen parameters - the gaugino mass parameters Mi, M2 and 
M3; the squark mass parameters for the first two generations Mq, Mu and Mr,, the third genera- 
tion squarks M3Q, Mbr, Mtr\ the corresponding slepton mass parameters M^, Mr, M^l, Mtr, 
the trilinear couplings At, and A^', and the higgs sector parameters Ma and tan /3, the ratio 
of the vacuum expectation values of the two Higgs doublets. 

The best reach in superparticle masses at the LHC is expected in the channel with two or more 
hard jets and missing energy |18[| which is the characteristic signature from gg and qg production. 
In particular, the simplest decays of the gluino and the squark, viz. g — > qqxi and q — )■ qx^ result 
in four and three-jet states with large missing energy for the two production processes respectively. 
Also, the sub-dominant qq* and qq production processes would result in two-jet final states. There- 
fore, looking for 2-4 jets with missing energy is known to be the best channel for SUSY searches. 
However, in the case that the first two generations of squarks are not accessible at the LHC (as can 

HQ 

indeed be motivated from suppression of flavour changing neutral currents )[|9|, 110|], the power of 
these searches would be dramatically reduced. Most importantly, the t-channel g^-type processes 
which contribute largest in the standard cMSSM analysis, would be severely suppressed because 
of the miniscule fraction of b and t-quarks in the proton. Similarly, other t-channel processes 
leading to qq, qq* are also suppressed. The dominant production processes are therefore gg and 
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s-channel qq* (q = ii,2, &i,2)- 

Since the limits on the mass of third generation squarks in a cMSSM-based analysis follow sim- 
ply from limits obtained from squark and gluino production processes, they cannot be considered 
truly indicative of the limit on stop and sbottom masses. We therefore reinterpret the ATLAS limits 
in the pMSSM model where the first two generations of squarks and all sleptons are decoupled. 

Considerable interest has also grown in recent times in SUGRA with non-universal high-scale 
masses. A high-scale parametrisation has the advantage that the masses of several particles are 
obtained naturally through renormalisation group (RG) running. Here too, we focus of situations 
where the third family sfermions are within the reach of the LHC, while the first two families are 
heavy. This, among other things, helps in a natural suppression of FCNC. The advantage of this 
scenario will be to allow us to investigate the effect of a low-mass slepton sector without requiring 
its full pMSSM parametrisation. 

We base our study on the data from the ATLAS experiment for signatures with jets and missing 
energy with and without b-tagged jets [|lU ll3l-ll6n. The ATLAS analysis assumes an mSUGRA- 
type unification for the interpretation of its data. As mentioned above, the limits from this analysis 
cannot be applied to the third generation squarks in a model independent way. The results have 
been interpreted in terms of a high-scale non-universal model in nllu . However, its dependence 
on mSUGRA based mass hierarchies (e.g. the lighter stop ii is always right-handed, the lightest 
neutralino is mostly bino-like etc.) hampers full understanding of the implication of the experi- 
mental data on the third generation squark sector. We therefore perform a more detailed study by 
performing a low- scale pMSSM analysis with a scan over the physical stop/sbottom masses. We 



also include the case of stop decay via the flavour- violating decay t — )• ex? lll8l - l20h when all other 
decays are forbidden by kinematics. 

The phenomenology of third generation squarks has also been studied in various scenarios [121 - 



3411 . It should be noted that while this manuscript was being prepared, similar studies on related 



issues [|35 



3611 have also appeared in literature. 



The paper is structured as follows: in the next section, we describe our simulation strategy 
and demonstrate the degree of agreement with the original ATLAS analysis. In section Unl we 
describe the phenomenological modelling of the stop and sbottom sectors and present the results 
corresponding to this scenario. The case of various high-scale non-universal scenarios is described 
in section HVl We present our conclusions in section IVl 
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II. CMSSM: SIMULATION OF THE SIGNAL AND ATLAS EXCLUSION CURVES 



We simulate the signal using Pythia 6.4 OTII and all strong production cross sections are 



normalised to their next-to-leading order (NLO) values as obtained from Prospino 2.1 11381] . The 
renormalisation and factorisation scales are set to the average of masses in the final state of the 
hard scattering process. We follow the detector acceptance region for all reconstructed objects and 
apply all the cuts as described in the ATLAS papers. 

The full set of identification and acceptance cuts is as follows: 

• Electrons: (1) px > 20 GeV (2) \r]\ < 2.47 (3) Sum of pr of particles within AR = 
a/ (A?7)2 + (A0)2 < 0.2 should be less than 10 GeV (4) Event vetoed if electron found in 
1.37 < \ri\ < 1.52. 

• Muons: (1) pr > 20 GeV (2) \r]\ < 2.4 (3) Sum of pr of charged tracks within AR = 



v/(A?7)2 + (A0)2 < 0.2 should be less than 1.8 GeV. 

Jets: (1) Formed using Anti-kt algorithm from Fastjet with parameter R = 0.4 (2) p^ > 20 
GeV. 

b-tagged jets: A jet is b-tagged if a b-hadron falls within a cone of radius R from a jet. We 



have c 



in im 



lecked that this reproduced the 50 % tagging efficiency for tt samples as mentioned 



• Missing transverse energy is calculated by summing over the pt of all objects and all stable 
visible particles not belonging to any reconstructed objects but falling within I77I < 4.9 with 
> 0.5 GeV. 

To account for detector effects, we smear the momenta of leptons and jets obtained from the 
Monte Carlo generator according to 

he values of a and b are (0.1 1, 0.007) for electrons, (0.03, 0.06) for muons and (1.0, 0.05) for jets 



[I52I1 . After smearing, we apply cuts used in each ATLAS analysis under consideration. 

A cross-check of our analysis is the correct reproduction of the missing transverse energy 
(MET) and effective mass (Mes) distributions and consequently the reproduction of the ATLAS 
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exclusion curves 11131 1 1411 in the context of mSUGRA. The online resources for the jets+MET anal- 
ysis for 35 pb~^ nun provide both the benchmark points used in the scan as well as the efficiencies 
and cross sections at each of the points. We use this information to verify the correctness of our 
simulation. We present in Fig.[Tl the final exclusion curves for the jets+MET analysis at 165 pb~^ 
which are obtained using the same acceptance and smearing parameters. We find that the "true" 
ATLAS curve lies in between our leading order (LO) and NLO curves in all cases. 

Since we aim to determine the limits on third generation squarks, the limits on b-tagged events 
are of particular importance and are expected to provide much stronger limits. Therefore, reliable 
modelling of b-tagging is of prime importance in looking at these signals. The work reporting 
the analysis of b-jet+MET describes the performance of the b-tagging algorithm as having an 
efficiency factor of 50% for a tt sample. We reproduce this number by this simple algorithm: 
we first form jets using the anti-Zcj- algorithm using Fastjet 2.4 ll39ll with the radius parameter 
i? = 0.4. A jet is assumed to be b-tagged if a b-hadron is found within a distance R from its axis. 
Since the correct reproduction of MET and M^s has been verified from the non-b-tagged samples, 
we can see that this algorithm gives a reasonably correct b-jet tagging by looking at the bottom- 
right panel in Fig.[TJ Here too, we find that the LO and NLO curves encapsulate the ATLAS curve 
reliably. 

In all the above cases, we find that the LO curve only slightly underestimates the ATLAS 
limits. In the worst case, the difference between the LO limits and the ATLAS curve is within 
20%. We therefore mostly present the LO mass limits in the subsequent study. Our LO results 
with simplistic detector simulation do not differ by more than 20% from what a full re-analysis 
of the data, including detector responses, would give. Our cross-checks on the mSUGRA results 
convince us such limits are adequate for putting across our main point, given the uncertainty of 
detector effects. 

We use the results on jets+MET for 1.04 fb"^ JlJ and b-jets+MET resuhs for 0.83 fb"^ Jm] 
for our analysis. As mentioned above, our simulation includes all the cuts in each of the channels 
under consideration. The benchmark points for our analyses are obtained using SUSPect 2.41 
l 4ol] . To obtain the exclusion curves, we use the value of cross section times acceptance provided 
by the ATLAS analysis. The values corresponding to different channels taken from these ATLAS 
analyses are summarised in Table HI The names of the signal regions are the same as those used in 
the respective ATLAS analyses. 

We concentrate here only on signals without leptons in the final state. The primary reason for 
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FIG. 1 : Comparison of our exclusion curves with the ATLAS exclusion curves 
sent: 2-jet with 165 pb~^ (top-left), 3-jet with 165 pb~^ (top-right), 4-jet with 165 pb~^ (bottom-left) and 
b-jet with 0-leptons with 35 pb~^ (bottom-right). The true ATLAS curve lies between our LO and NLO 
contours in all cases. 



this is that we wish to investigate the third generation squark sector in as model independent a way 
as possible. Leptonic states generally result from decays on gauginos into gauge bosons or into 
sleptons which then decay into leptons. For a pMSSM study based on leptonic signatures, it would 
therefore be imperative to also include a completely general gaugino sector as well as a low-mass 
slepton sector. Since adding a completely phenomenological slepton sector means adding five new 
parameters which complicate the analysis beyond too much, our low-scale analysis deals with a 
decoupled sleptonic sector. It is, in principle, possible that the limits obtained in the decoupled 
slepton limit are diluted when decays into sleptons (and hence leptons) are possible. We include 
the possibility of a low-mass slepton sector in section|IVl where we use high-scale parametrisation 
and allow the RG running to determine the masses in the slepton sector. However, we shall see that 
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allowing low-mass sleptons do not make significant difference to the limits for signatures based 
on jets and missing energy. 

We perform our analysis retaining the cMSSM-like gaugino mass pattern Mi : M2 : M3 = 1 : 
2 : 6. In general, for most gaugino mass patterns, we expect our results to remain fairly unchanged 
since our signal does not depend strongly on particles obtained from intermediate decays in SUSY 
cascades. However, we explicitly comment at the end of the paper on the extreme cases of gaugino 
mass patterns that would be likely yield results drastically different from ours. 

It is also possible to ask what fraction of the pMSSM phase space is ruled out by current data. 



For the effect of the experimental limits on the full pMSSM space, we refer the reader to 114 iL 1421] 



The effects on the cMSSM parameter space are addressed in, for example^^ 



interpretations of the recent LHC data on SUSY have been discussed in ||4| 



M42t4 
46^491]. 



2|; 45|] whereas other 



Channel 



a X acc (fb) 



2 jets + MET 


24 


3 jets + MET 


30 


4 jets + MET (Meg =1 TeV) 


32 


1 btag + Mcff > 500 (3JA) 


288 


1 btag + Mcff > 700 (3JB) 


61 


2 btag + Afcff > 500 (3JC) 


78 


2 btag + Mcff > 700 (3JD) 


17 



TABLE I: The values of cross section times acceptance from ATLAS analysis used for 
limits. The first set uses 1.04 fb~^ of data 1 15] whereas the second uses 0.833 fb~^ of data 



applying exclusion 
116]. 



III. PARAMETERISATION OF THE STOP AND SBOTTOM SECTOR 

We work in the plMSSlVI framework where the parameters are assigned at the low-scale. The 
program SuSpect is used to ensure that electroweak symmetry breaking has correctly taken place 
and the spectrum is consistent. To start with, we retain the clVISSlVl-like gaugino mass ratios, 
correspond to Mi : M2 : M3 ~ 1 : 2 : 6 among the U(l), SU(2) and SU(3) gaugino masses. 
The squark masses of the first two generations and all the slepton masses are set to 2 TeV which 
is beyond the reach of the 7 TeV LHC run. The stop and sbottom sector can each be described 
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by three parameters - the two mass eigenstates and the mixing angle. We use these as the input 
parameters for the scan and use the diagonalisation to determine the left and right handed mass 
parameters of the pMSSM. The stop sector requires three parameters — the masses M^^ and M^^ 
and the stop mixing angle 6^. Using the stop mass-squared matrix diagonalisation condition 

M? \ I A/2 mtXt\ , I cosOi sinoA 

Ml J \mtXt j \-sin^4 cosOi j 

where Xt = At — ji cot /3, we can use the low-scale masses Mi^ , M^^ and 6^ as the input parameters 
which uniquely determine At given fj, and tan /3. The left handed sbottom mass is expected to be 
close to the left handed stop mass since they are derived from the same parameter (M^q) in the 
low-scale pMSSM model. The right handed sbottom mass and Aj, can then be set depending on 
the requirement of M^^, M^^ and O^. 

Assuming that the third generation squarks and the gluino are the only strongly charged super- 
particles accessible at the LHC, we investigate in particular, the following cases: 

• Case A: sinOt = 0.99 i.e. ii ^ iji is the lightest squark. This is commonly the case in 
cMSSM models. We also set A/^^ ~ A% ~ Af^^ = Mf^ + 500 GeV, which makes the 
sbottom sector somewhat heavier than the lighter stop. We scan the parameter space in 

M = Mr . 

• Case B: sm% = 0.01 i.e. ii ~ tL, ~ 6i ~ 6^. The lightest stop and sbottom are nearly 
degenerate and mostly left handed. Mp^ ^ M^^ = Mi^ + 500 GeV. We scan the parameter 
space in M = M^^ = M^^. 

• Case C: 6i ~ Ir is the lightest squark. Af^^ ~ Mt^ ~ M^^ = M^^ + 500 GeV with 
sin 6t = 0.70 and sin 6^ = 0.99. We scan the parameter space in Af = Af^^. 

• Case D: Af^-^ ~ Af^^ ^ Af^^ — ^ba- This is the case of maximal mixing in both stop and 
sbottom sectors. The most stringent limits on the light third family scenario will arise for 
this particular case, as it allows all four squarks to be produced with similar cross sections. 
The scan in this case is over the common mass of all the third generation squarks. 

In each case, we are now able to perform a scan over the Af — M2 plane, where Af is the mass 
of ii in cases A, B and D and bi in the case C. The Tevatron reach for searches in the ii — )■ ex? rule 



out stop masses up to 180 GeV bOll . We therefore start our search at Af^^ = 200 GeV and scan up 
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Parameter Scan range 



100 - 2000 GeV 
M2 150 - 600 GeV 

tan /3 5, 10, 40 

fj, -200, 200, 500, 1000 GeV 



TABLE II: The parameters of the scan in M^_^-M2 space. 

to 2 TeV. Also, since we assume gaugino mass unification, we use the chargino mass limits from 
Tevatron (M-+ >164 GeV) to start our scan at M2 = 150 GeV and vary M2 up to a value of 600 
GeV which would correspond to gluino mass of 1.8 TeV and therefore cover the entire range of 
masses reachable at the 7 TeV run of the LHC. We fix Ma = 400 GeV and perform this scan for 
12 combinations of 3 values of tan /? and 4 values of /i, which are listed in Table UIl 

The hierarchy among Mi, M2 and /i determines the composition of the neutralino and chargino 
sector and therefore has strong effects on the limits. In particular, for cases with low values of fi, 
the lightest neutralino has a significant higgsino fraction. When M2 becomes large enough, the 
states Xi and xt become higgsino dominated and their masses remain close to the value of //. In 
this case, the masses and compositions of Xi^ Xi X2 affected by changes in the value of 

M2, so long it is considerably larger than fi. For large fi, on the other hand, the allowed decays 
of the squark will depend strongly on M2 up to very large values. Therefore, we expect that for 
low squark masses and fi low with respect to M2, the exclusion contour is relatively insensitive to 
M2. This can clearly be seen in figure [2] and best illustrated in the fourth panel corresponding to 
Case-D. Here, the production cross section is high because all four third generation squarks are 
degenerate. Moreover, it can clearly be seen that /i = -200 and 200 both lead to large exclusion 
in M2 for small masses. The effect is similar also in the panel corresponding to Case-B. The third 
panel, corresponding to sbottom being the lightest shows minimal change with changing fi. This 
is mostly because the decay bi — )■ bxi is always open irrespective of fi due to the small mass of 
the b-quark. Thus, the sensitivity to /i is reduced. 



The channel with two b-tagged jets and high missing energy (called 3JD in nlw ) results in 
most stringent limits in all cases. The cross section ruled out in this case (0.017 pb) is similar to 
the the number for 2 jets + MET (0.024 pb), however, the latter does not contribute significantly 
to the parameter space ruled out. The reason for this is two-fold: first, the effective mass cut of 1 
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FIG. 2: The dependence of the exclusion curves on different values of f-i for tan /3 =10. The least depen- 
dence is for the case where the sbottom is the lighter squark. The shaded region at the top in the fourth 
panel corresponds to stop LSP and is therefore ruled out. The x-axis refers to M^^ for Cases A, B and D 



and Mr for Case C. 

01 



TeV for the jets+MET final states has a much smaller efficiency for small stop masses. Secondly, 
the missing energy cut for 3JD is a smaller, fixed number (130 GeV) as compared to the M^s 
dependent one (MET > O.SMoff) in the 2 jets + MET case. The reduction of effeciency when 
demanding 2 b-tagged jets is not very severe in gg production events due to four b-quarks in the 
final state. Even in the regions where exclusion is dominated by gg*-type production, we still find 
that the 3JD channel results in the best exclusions. 

To estimate the effect of enhancement due to NLO corrections in the production cross section, 
we present the comparison of LO and NLO curves for each case with tan (3 = 10 and fi = 500 GeV 
in Fig.[3l The large k-factors (~ 2.5) in most of the parameter space result in much stronger limits 
from the NLO curves. However, taking note of the results of our cMSSM limits in Section [III 
where the LO limits are closer to actual ATLAS limits, we take the conservative approach of 
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FIG. 3: Comparison of LO and NLO exclusion curves for Cases A, B, C and D, illustrated for values 
tan (3 = 10 and fi = 500 GeV. The x-axis refers to mass of the lightest third generation quark in each case, 
as discussed in the text. 

presenting LO limits for our study. 



Results: Case A 

This is the case closest to mSUGRA-type models where ii is the lightest squark. The primary 
production processes in this case are gg and tit*. Since our scan starts with M2 = 150 GeV, i.e. a 
gluino mass of 450 GeV, the decay g — )■ ttxi is always open and forms the dominant decay mode. 
In the low mass regions which are probed by the LHC data, the dominant decay mode of the 
ii depends on the mass hierarchy of ii, and xf. Since we are working in R-parity conserving 
models, we disallow the region where the lighter stop is the lightest supersymmetric particle (LSP) 
i.e. the region where M^^ < M^o. Following this, the hierarchy M^o < M^^ < M-+ results in 
the case of stop NLSP (next-to-LSP). The dominant mode is ii tx\ if M^^ > Mt + M^o. 
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FIG. 4: Exclusion curves for Case-A (stop lightest scenario). The x-axis refers to the mass of ti squark. 

Otherwise, we expect ti to decay via three-or four-body decays or via the mode ii cx^. For 
this work, we assume that this last mode dominates over the three or four body decays. Finally, in 
the case where M^^ > Mb + M^+, the decay into bxt is also open. For larger stop masses, decays 
into other neutralinos, second chargino or gluino may also open. 

The effect of various values of tan f3 for two values of are shown in Fig. IH The exclusion 
curves show minor dependence on tan /?. We find that for low values of /i, low tan /3 results in a 
larger reach in M2 whereas high tan (3 results in a larger reach in M^. This trend is reversed for 
large values of /i, as can be seen in the panel corresponding to = 500 GeV of Fig.|4l However, 
it must be reiterated, that this variation cannot be considered experimentally significant due to the 
uncertainties on the exclusion curves in our analysis. 



Results: CaseB 

In this case, the ii and bi form a degenerate pair of lightest squarks. They are both primarily left 
handed and therefore have an enhanced coupling to wino-like states. Again, due to the requirement 
of neutral dark matter candidate, we disallow any region with stop or sbottom LSP. The decays of 
the lighter stop are similar to those in case A. The decay of the lighter sbottom into bxi is almost 
always allowed and will form the dominant decay for most of the low mass region. In cases of 
large fi, where M^o ~ M2, the decay 61 — )■ will dominate over 61 — )■ bxi similarly for stop 
decays. The decay bi txt is relatively disfavoured due to large top mass. The gluino decays 
dominantly to bbxi in the region Mg < M^^ and to bbi otherwise. The decays to corresponding 
top-sector are again disfavoured due to large top mass. 
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FIG. 5: Exclusion curves for Case-B (left-handed degenerate stop and sbottom scenario). The x-axis refers 
to the common mass of bi and ii squaife. The shaded region at the top in the fourth panel corresponds to 
stop/sbottom LSP and is therefore ruled out. 
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FIG. 6: Exclusion curves for Case-C (sbottom lightest scenario). The x-axis refers to the mass of bi squark. 



The dependence of the exclusion curves on tan (3 is shown in Fig. [5l For the case of ji = —200 
GeV and tan /? = 5, we find that ti and hi up to 300 GeV are ruled out for gluino masses up to 
700 GeV. The region just below the ii or 61-LSP region is still allowed, as the near-degeneracy of 
their mass and the mass of the LSP results in a low missing energy and Meg spectrum which does 
not satisfy the hardness cuts imposed. 



Results: Case C 



This case considers the situation where both the stop states are heavier than the lightest sbottom 
state. The primary production processes are gg and hih\. As before, we disallow the region 
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FIG. 7: Exclusion curves for Case-D (maximal mixing scenario). The x-axis refers to the common mass of 
all the third generation squarks. The shaded region at the top corresponds to stop LSP and is therefore ruled 
out. 

< M^o. The region M^o < M-^^ < M-+ corresponds to a 61-NLSP with the dominant 
decay bi — )• bxi ■ The gluino dominantly decays via bbxi in the region Mg < M^^ and to bbi 
otherwise. This case is the closest to the scenarios considered by the ATLAS collaboration for 
the interpretation of their b-jet and missing energy searches. They split their analysis into a case 
where they disallow any three body decays of the gluino via g — t- bbxi ^^^^ of ^ simplified 

model where there are no two-body decays but all decays are via this channel. According to the 



first case, they rule out gluino masses up to 720 GeV for 61 masses up to 600 GeV 111611 . As can 
be seen from Fig.[6l for bi = 600 GeV, we rule out Mg < 570-660 GeV for = -200 and Mg < 
600-660 GeV for /x = 500 GeV. 



Results: Case D 



In the maximal mixing scenario in both the stop and sbottom sector, all four squarks of the third 
generation have nearly degenerate masses. Therefore the production cross section is maximum for 
this scenario and the limits are strongest. The decay scheme for the sbottom is same as in Case 
B whereas for the stop, it is the same as in Case A. The gluino can now decay both via stops 
or sbottoms, but the large mass of the top means it decays preferentially via bbi channels. As 
expected, low fi results in large exclusion in M2 for low squark masses. The dependence of the 
limits on tan (3 is shown in Fig. |7J The case for fi = —200 GeV (and similarly /i = 200 GeV) 
result in an exclusion of third generation squark masses of 280 GeV for all allowed values of M2. 
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This case is similar to the one considered in [13511 . where, for naturalness requirements they require 
iji, II and to be degenerate and assume a Higgsino LSP. Their limits on the mass of the third 
generation squarks lie between 200-300 GeV. 

The exclusion limits for the case of /i = 500 GeV does not show any exclusions that are 
independent of gaugino (and hence gluino) masses. The third-generation squark masses are com- 
pletely un-constrained for Mg > 1 TeV. Whereas, the approximate requirement for naturalness 
that < 500 GeV would translate into Mg > 600 GeV 

In conclusion, among the four cases discussed above, only the case with degenerate third- 
generation squarks and low-/i leads to mass limits independent of gluino mass - that of 280 GeV. 
In most other cases, we find that limits depend strongly on the composition of the neutralinos 
and charginos. The case where only ii is accessible is the least constraining, mainly due to low 
production cross sections compared to the other cases. For the case where the LSP is a almost 
pure Bino (high-//), M^^ = 200 GeV is ruled out for a gluino mass less than 570 GeV. Taking into 
account all values of /i and tan /3, Mi^ = 200 GeV is ruled out for gluino masses in the range 
570 - 720 GeV. 

The case of lightest third generation squark being bi is the most insensitive to variations of both 
fi and tan (3. For this case, our limits are consistent with ATLAS's own interpretations within 10%. 
M^^ = 200 GeV is ruled out for gluino masses between 680—820 GeV. The case of degenerate left- 
handed squarks rules out M^^ = M^^ = 200 GeV for gluino masses in the range 900 — 1050 GeV. 
And finally, the case with all squarks degenerate rules out j = 200 for Mg < 900 GeV in the 
worst case, and for all gluino masses in the best case. 



IV. HIGH-SCALE NON-UNIVERSAL SCALAR SCENARIOS 

Besides the low-scale study done in the previous section, it is also possible to perform a scan 
over high-scale parameters. The advantage of a high-scale analysis is that the hierarchy of the 
particles is uniquely and consistently determined from the renormalisation group (RG) running of 
masses to low scale from given high-scale parameters. We use the simplification afforded by this 
model to include the effects of sleptonic sector in our analysis. 

It is possible that the limits described in the previous section are diluted if slepton masses are 
allowed to be light. This is because the gauginos would then decay predominantly to sleptons 
resulting in leptonic final states which would be discarded since all the analyses considered here 
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have a lepton veto. Including the full lepton sector in the low-scale pMSSM requires five more 
parameters and makes a general study far more complicated. We therefore leave a fully model 
independent investigation of interpreting the ATLAS limits involving a low mass slepton sector to 
a future work. However, we partially answer the question as to whether the limits are diluted by 
studying some illustrative cases, as described below. 

Even though the soft-scalar masses may in principle take separate values, the constrains from 
flavour changing neutral currents (FCNC) from meson decays dictate that the first two generation 
squarks remain degenerate. Similarly, absence of decays like /i — > 76 means that the masses of 
first two generations of sleptons also have to degenerate. Therefore we can consider three schemes 
of non-universality: 

• Case HA: Third generation squarks are lighter than all other sfermions. This will lead to a 
hierarchy similar to Case A in the preceding analysis. 

• Case HB: Third generation squarks and sleptons are lighter than all other sfermions. This 
leads to light staus and tau-sneutrinos. Possibly, this would lead to br final states which have 



been studied in 115 ih . 



• Case HC: Third generation squarks and all sleptons are light. Comparing the limits in 
this case to those in Case HA will answer the question of dilution of limits due to leptonic 
signatures. 

As in the previous analysis, we retain the cMSSM-like gaugino masses. We now have two 
mass scales in the scalar sector - the scale of heavy, decoupled particles (Mheavy = 2 TeV) and the 
scale of light scalars (Mught). We once again assume that the gaugino sector follows the universal 
structure and we set the higgs mass parameters for the two higgs doublets to be same as Mught. 
The parameters are all set at the grand unified theory (GUT) scale and the TeV-scale values are 
determined by RG running using the program SUSPect. The exclusion curves obtained for the 
three cases are shown in Fig. [8] for a value of tan /3 = 10. The results for values of tan/3 = 
5, 40 are similar. We do not see any significant difference among the exclusion curves for the 
various high-scale cases. This can be interpreted as the robustness of the zero-lepton signals 
against different slepton masses and justifies the assumption of decoupled slepton masses made in 
the previous section. 
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FIG. 8: Comparison of exclusion curves from three high-scale non-universal scenarios. All cases agree with 
each other within statistical uncertainties. 

V. CONCLUSIONS 

We have investigated the consequences of the recent ATLAS data in channels with (b-)jets and 
missing energy on the limits on the mass of the third generation squarks. We work in the pMSSM 
framework, with TeV-scale parameters, without requiring a high-scale breaking scheme. For ob- 
taining relatively model independent limits on the third generation squark masses, we decouple 
the first two squark generations as well as all sleptons. We also explicitly show that decoupling of 
sleptons is not likely to affect the limits as long as we work with 0-lepton signatures. We find that a 
stop of mass 200 GeV can be ruled out for a gluino mass of 570 GeV in the least constraining case 
whereas a stop of mass 500 GeV is allowed for gluino masses upward of 450-880 GeV depending 
on the structure of the third-generation squark sector, and the parameters /i and tan f3. In the case 
where all third generation squarks are degenerate, we can rule out masses less than 280 GeV for 
< 200 GeV, independent of the gluino mass. 

One may also question the assumption of the gaugino mass pattern of Mi : M2 : M3 ~ 1 : 2 : 6. 
We can expect the case of right-handed sbottom being the lightest of the third generation to be 
fairly independent of the assumption of gaugino mass pattern since the only dominant decays are 
bi bxi and g — )■ bbjcf, and both are always allowed (except for very compressed spectra). We 
would expect significant deviations from the stop limits when, for example, the decay ii — t- bxt 
is largely inaccessible because both M2 and /i are so high that the charginos are generally heavier 
than the stop. In this case, the dominant decay for most of the parameter space would be ii — txi. 
In the current study, regions where the txi decay was kinematically disallowed was still largely 
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covered by the bxf decay, thus leaving only a small region of parameter space corresponding to 
the flavour violating decay. However, in the absence of the decay into a chargino, one would 
need to examine in detail, the relative strengths of the (highly model-dependent) flavour-violating 
decay ii — cxI and the three-body decay of the stop. In other cases, where the stop still decays 
via standard channels, we do not expect significant deviations from our limits. 
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